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7 The charging of spinning insulating objects by plasma and photoemission 

8 is studied with the particle-in-cell method. Unidirectional photon flux, dif- 

9 ferent angular velocities of the object, and different plasma flow speeds are 

10 considered. Photoemission can lead to a positive total charge and electric dipole 

11 moment on the object. The spinning of the object redistributes the surface 

12 charge. The total object charge oscillates in time with the period matching 

13 the period of the object rotation. Plasma potential and density in the vicin- 

14 ity of the object oscillate with the same frequency. The plasma is rarefied 

15 close to a positively charged object, and the density wake diminishes when 

16 the object is charged negatively. The time averaged charge depends on the 

17 angular velocity of the object. 
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1. Introduction 

18 The charging of an object by plasma is one of the basic problems in space and plasma 

19 physics. If only plasma currents are considered, the charge on the object is usually negative 

20 [Svenes and TY0im, 1994; Vladimirov et al, 2005; Miloch et al, 2007]. The plasma flow 

21 introduces asymmetry to the object's charging and gives rise to wakes in the plasma den- 

22 sity and potential [Vladimirov and Nambu, 1995; Melands0 and Goree, 1995; Vladimirov 

23 and Ishihara, 1996; Ishihara and Vladimirov, 1997]. This asymmetry is more pronounced 

24 for insulating than for conducting objects [Miloch et al, 2007, 2008a]. Photoelectric cur- 

25 rent can lead to positively charged objects [Shukla and Mamun, 2002; Vladimirov et al, 

26 2005] . If photoemission is due to a directed photon flux, the electric dipole moment can 

27 develop on insulating objects. There is a significant difference between the wake of the 

28 negatively and positively charged object, with a strong density rarefaction for the latter 

29 [Miloch et al, 2008b]. 

30 The understanding of the charging of an object in a complex environment with sunlight 

31 and plasma flow is of concern for the operation of spacecrafts or sounding rockets [Svenes 

32 and Tr0im, 1994; Roussel and Berthelier, 2004]. The charging of insulating components 

33 of such objects can result in strong potential differences between the shadow and sunlit 

34 sides. Changes of the plasma parameters in the vicinity of the object need to be accounted 

35 for when analyzing the instrument data [Lai et al., 1986]. The problem can be further 

36 complicated by the spinning of an object, whether it is an intrinsic rotation of an asteroid 

37 or dust grain, or imposed for the attitude and stability control of a spacecraft or rocket 

38 [Lee et al., 2001; Kurihara et al, 2006]. 
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39 Analytical models for satellites spinning in sunlight demonstrated the development of 

40 potential barriers that decelerate photoelectrons and allow the sunlit side of the spacecraft 

41 to be negatively charged [Tautz and Lai, 2006, 2007]. These models assume the satellite 

42 to be much smaller than the Debye length, and replace the Poisson equation with the 

43 Laplace equation for vacuum. Moreover, they do not account for the plasma flow. If the 

44 plasma flow and intermediate sizes of the satellite with respect to the Debye length are 

45 considered, the problem becomes highly nonlinear. A self-consistent analytical model for 

46 such a problem is difficult to develop, and the numerical analysis seems appropriate. 

47 In this letter we discuss results from the particle-in-cell (PIC) simulations of spinning 

48 insulating objects in flowing plasmas exposed to unidirectional photon flux. The analysis 

49 is relevant for such objects as satellites, rockets or asteroids in space, boulders on lunar 

50 surface, or dust grains in experimental devices [Hordnyi, 1996; Fortov et al, 1998; Khrapak 

51 et al, 1999]. 

2. Numerical code 

52 The analysis is carried out in two spatial dimensions in Cartesian coordinated using the 

53 PIC numerical code described in detail by Miloch et al. [2007, 2008a, b, c]. The electrons 

54 and ions are treated as individual particles, with the ion to electron mass ratio being 

55 mj/mg = 120. The plasma density is n = 10^°m~^. The coUisionless plasma flows in the 

56 positive X direction, and three plasma drift velocities are considered = {0, 0.5, LSjC^, 

57 with Cs being the speed of sound for adiabatic ions and isothermal electrons: Cg — 



K.{Te + 'jTi)/M, where 7 = 5/3. The electron to ion temperature ratios are C = T^/Ti — 




59 



{5, 100}, where = 0.18 eV. 
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60 A circular object of radius r = 0.375 in units of the electron Debye length A^e is placed 

61 within a simulation box of size of 50 x 50 Xoe- Such an object can be understood as 

62 an intersection of a cylindrical object in a three dimensional system, and it is initially 

63 charged only by the collection of electrons and ions. To represent a perfect insulator, 

64 each plasma particle that hits the surface remains at this position for all later times 

65 contributing to the local charge density. The object spins throughout the whole simulation 

66 with angular velocity cu of 0.57r, 27r or alternatively Stt in units of rad/rj, where Tj is the 

67 ion plasma period. A directed photon flux is switched on at 39Tj, when it can be assumed 

68 that the surface charge distribution has reached a stationary level. The code is run for 

69 approximately 50Tj. 

70 The code allows for an arbitrary angle of the photon incidence on the object. When 

71 a photon hits the object surface, a photoelectron of energy 1 eV is produced at distance 

72 I — svAt from the surface, where s is an uniform random number s e (0, 1], A.t is the 

73 computational time step, and v is the photoelectron speed. Photoelectron velocity vectors 

74 are uniformly distributed over an angle of tt and directed away from the surface which is 

75 in accordance with Lambert's law. In this study the incoming photons are usually aligned 

76 with the direction of the plasma flow, i.e., a = 0°. As a control case, we also consider 

77 the angle of 90° between the incoming ions and the plasma flow {a — 90°). The photon 

78 flux is = 1.25 X lO^^m^^s^^. The scheme of a typical numerical arrangement is shown 

79 in Fig. 1. Points close to the object's surface, that are labeled with numbers, indicate 

80 probes for potential variations that are described in more detail in Results section. 
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3. Results 

81 With the onset of radiation, the total charge on a spinning object becomes more positive 

82 and starts to oscillate in time. The period of oscillations matches the period of the full 

83 rotation of the object, see Fig. 2. The charge variations are large for slowly spinning 

84 objects. The mean charge depends on the plasma flow velocity. It reaches lower values for 

85 faster plasma flows Vd, and is more positive with higher angular velocities of the object cu. 

86 The mean charge values for different Vd and cu are summarized in Table 1. Since there is 

87 little difference in the charging characteristic for ^ = 5 and ( — 100, only the results for 

88 ( — 100 are shown in Fig. 2 and Table 1. 

89 Photoemission due to unidirectional photons leads to the development of an electric 

90 dipole moment on the object. The electric dipole moment is initially antiparallel to the 

91 photon direction, and it co-rotates with the object. It vanishes at the certain angle between 

92 the incoming photon direction and the dipole moment, which increases with increasing cu. 

93 The electric dipole moment antiparallel to the photon direction reappears after the full 

94 rotation of the object. 

95 In flowing plasmas, the wake in the plasma density forms behind a spinning object. This 

96 wake oscillates in time. The region of rarefied plasma density has large spatial extent when 

97 the object is positively charged, while for a negatively charged object, the plasma density 

98 in the wake can be enhanced and the ion focusing observed. The oscillatory nature of the 

99 wake is demonstrated in the ion density plots in Fig. 3a) . The transition between the ion 

100 wake and the enhanced ion density is asymmetric. The edge of the wake is distorted by 
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101 the local enhancement in the ion density. This distortion grows in time and reduces the 

102 wake size. 

103 Similar rarefication in the ion density and scenario for the wake closure are observed 

104 for stationary plasmas, with the difference that the ion density is rarefied close to the 

105 positively charged side of the object and the wake does not form. 

106 Without photoemission, a symmetrical sheath forms around the object, while for flowing 

107 plasmas the wake behind the object is observed. The potential distribution in the vicinity 

108 of the object is governed by the photoemission and associated electric dipole moment. The 

109 oscillations of the total charge alter the potential distribution, and for flowing plasmas, 
no the wake in the potential is highly distorted, see Fig. 3b). Since many instruments on 

111 the rocket and satellite payloads are installed on extended booms, it is vital to examine 

112 to what extend a spinning insulating body modifies the surrounding plasma. In Fig. 4 we 

113 show potential variations at four point- like probes at distances d — Xoe from the object 

114 surface in the direction parallel and perpendicular to the photon and plasma flows. The 

115 orientation of each of the probe relative to the plasma flow is shown in Fig. 1. For all 

116 probes, we observe oscillations with the periods equal to the rotation period of the object. 

117 The amplitudes of oscillations on the side of the object that is charged predominantly 

118 positive (i.e., probes 1 and 2) are moderate as compared to probes 3 and 4. For the latter 

119 probes, other strong, repeatable components are present in the signal. The amplitude of 

120 potential variations for probes 3 and 4 are twice as large as for other probes. 

121 For temperature ratio C = 5, the total charge on the object is less negative than for 

122 C = 100 due to the ion mobility. The charging characteristic with photoemission is similar 
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123 for both cases, with the negative part of oscillations being more negative for colder ions. 

124 The ion wake is smaller and less pronounced for warmer ions, but the analyses of potential 

125 oscillations close to the object are similar for both temperature ratios. 

126 For the control case of the photons incidence angle of 90° with respect to the plasma 

127 flow, the wake is spatially different from the case for a — 0°. However, the principal 

128 mechanism and oscillations of the total charge in time are the same. The mean charge is 

129 slightly higher for a — 90° than for a = 0°. 

4. Discussion and Conclusions 

130 With the onset of the photoemission due to unidirectional photon flux, an electric 

131 dipole moment develops on the object. Most of the positive charge is localized on the 

132 illuminated side of the object, while the opposite side is predominantly negatively charged. 

133 For a spinning object, the photoelectric current neutralizes negatively charged regions, 

134 but the photoemission rate can often be too low to charge them positively. Only after 

135 approximately one rotation period of the object, photoemission can lead to the positive 

136 total charge on the object and recovery of the electric dipole moment. This suggests 

137 that the initial charging of the object by photoemission is robust, and that the charge 

138 redistribution on the object's surface is insufficient to compensate for the residual charge 

139 with the spinning rates considered in this study. 

140 This picture should be complemented by the ion and electron dynamics close to the 

141 object. In case of the plasma flow, the wake is formed. For spinning objects, this wake 

142 is distorted on one side. Spinning, negatively charged side of the object accelerates ions 

143 towards the wake and enhances the plasma density locally. The wake is further distorted 
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144 due to the enhancement in the plasma density and neutrahzation of the positively charged 

145 regions on the surface. 

146 A spinning insulating object excites waves in the system with the wave frequency match- 

147 ing the frequency of the object rotation. While these waves are damped further away from 

148 the object, they are conspicuous close to the object's surface. Since the plasma density is 

149 inhomogeneous and the wake oscillates, the wave propagation is complicated. On the side 

150 predominantly positively charged, the variations in the potential distribution are smooth 

151 with a single dominant frequency observed. On the negatively charged side, where the 

152 wake oscillations modify the plasma to a high degree, the potential variations are larger, 

153 and modulations of the main wave are observed. 

154 Previous analytical works by Tautz and Lai [2006, 2007] did not include plasma dynam- 

155 ics around the object, and hence stationary solutions could be obtained. In the limiting 

156 case of a small object with respect to Xoe, one can expect such results with the conspicu- 

157 ous potential barrier for photoelectrons, especially for fast spinning insulating objects. We 

158 have shown, however, that in the regions where the plasma dynamics around the object is 

159 important, the result will be different, and the photoemission and spinning of the object 

160 break the symmetry of the object charging even without the plasma flow. In particular, 

161 the satellite or rocket instrument readings can be influenced by the rotation, which has to 

162 be accounted for in data analysis. For a spinning satellite traveling through regions with 

163 different plasma densities and temperatures it will often be necessary to employ these two 

164 different models. 
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Figure 1. Scheme of a typical numerical arrangement. Points labeled with numbers corre- 
spond to the probes for potential variations in the vicinity of the object. The object rotates 
anticlockwise. 
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Figure 2. The charging of spinning insulating object during photoemission for different angular 
velocities to and speeds of the plasma flow Vd for C, — 100. The results are smoothed with a moving 
box filter for presentation. 
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Figure 3. Ion density (a) and electric potential (b) for a spinning insulating object during 
photoemission with to = 2tit^^ ^ Vd = 1.5Cs, and ( = 100. Different time instances are shown for 
which the total charge on the object is positive (top), the erosion of the wake begins (middle), 
and the total charge on the object is negative (bottom). 
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Figure 4. Potential variations for different probes at distances d — Xm the surface of the 
spinning object for Vd = l-5Cs, ou = 2tit^^, and ( = 100. The orientation of the probes with 
respect to the plasma flow and direction of photons is shown in Fig. 1. 



Table 1. The mean charge value on the spinning insulating object in the presence of photo- 
emission for ^ = 100. The charge is presented in units of elementary two-dimensional charge: 
qo — ey/n, where e is the elementary charge in a three dimensional system and n is a plasma 
density in a two dimensional system. 
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